Separating parallel conduction from two-dimensional magnetotransport in high 
mobility InP/InGaAs MOCVD-grown heterostructures 
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In this Letter, four-point magnetotransport of high mobiUty InGaAs/InP heterointerfaces is mea- 
sured from 1.6 K to 300 K and from to 15 T, and an analysis is shown whereby the mobiUty 
and density of the two-dimensional (2D) accumulation layer can be separately characterized from 
that of the parallel conducting dopant layer over all but a small intermediate temperature range. 
Standard magnetotransport regimes are defined as the temperature increases from 1.6 K to 300 K, 
namely quantum Hall (QH), Shubnikov de Haas (SdH), and Drude regimes (D), and in the QH 
and D regimes different analyses are applied to deduce densities and mobilities of both layers sep- 
arately. Quantitative conditions for the intermediate SdH regime are defined, within which both 
QH and D analyses fail. The density and activation energy of unintentional donors at the InP 
epilayer/substrate interface is deduced. At base temperature, QH minima are resolved down to 
B — 0.4 T at V — 20, revealing a mobility of /i = 160, OOOcm'^/Vs. The 2D system maintains this 
high mobility up to at least 40 K in this high quality structure. 

PACS numbers: 73.40.kp, 73.43.Fj, 73.43.Qt 



High efficiency quantum cascade lasers require high 
quality MOCVD InP/InGaAs heterostructure growth, 
which can be characterized by examining electrical trans- 
port at a single InP/InGaAs interface. Such character- 
izations of high mobility two-dimensional electron gas 
(2DEG) heterostructures are often obscured by paral- 
lel conduction arising from intentional or unintentional 
doping from an undepleted bulki"— , a second occu- 
pied subband^, an overdoped modulation doping layer—, 
or unintentional interface defects and impurities which 
can act as dopants in a localized layer. Various par- 
allel conduction models have been used to character- 
ize parallel conducting channels in both quantizedS and 
non-quantized^i^ magnetotransport of high mobility sys- 
tems. But to date, experimental characterizations of 
InP/InGaAs heterostructures are lacking, typically be- 
cause nominally undoped bulk layers nonetheless have 
moderately high background doping levels, and because 
interface defects and impurities induce localized donor 
states at the epilayer interface of unknown density. To 
date, parallel conduction characterizations have been 
performed principally at a single temperature, and tem- 
perature dependent characterizations for parallel chan- 
nels over a wide temperature range was rarely reported. 

In this paper, we present an analysis of magnetotrans- 
port data for an InP/InGaAs heterostructure from 1.6 
K to 300 to separately extract the carrier density and 
mobility of the parallel conducting Drude layer from the 
high mobility layer. Three magnetotransport regimes are 
labelled according to the conductivity regime of the high 
mobility layer, namely quantum Hall (QH), Shubnikov 
de Haas (SdH) and Drude (D). In the high temperature 
D range, a new technique is presented whereby the mag- 
netic field dependence of the conductivity over a large 
field range permits separation of the two parallel con- 
duction parameters. In low temperature QH range, the 
conductivity in the resistance minima reveals the conduc- 
tion parameters of the parallel layer. The high quality 
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FIG. 1: (a) Schematic configuration of the InP/InGaAs 
heterojunction. (b) Conduction band structure at the 
InP/InGaAs interface. Ep, Ad, and Ao are the fermi level, 
binding energy of interfacial dopants, and the first electron 
subband in the quantum well respectively. The cross stands 
for ionized donors and the wavefunction is corresponding to 
the first subband. 



of the heterostructure growth allows us to identify epi- 
layer interface states as the principle source of the parallel 
channel in this InP/InGaAs system, and a quantitative 
model for identifying the heterolayer conduction band 
offset from the data is discussed. 

The sample is grown by metalorganic chemical vapor 
depositionZ.iS and the structure is in Fig.[lja). The com- 
position of InGaAs layer is Ino.532Gao.467As. Neither the 
InP nor InGaAs layer was intentionally doped. Because 
of the high growth quality, it is known that the back- 
ground doping of the bulk layers is at most lO^^cm"'^. 
The interface between InP and semi-insulating InP sub- 
strate is known to be a possible additional source of 
donors for creating a parallel conducting channel, which 
is expected to dominate parallel conduction in such high 
growth quality structures. Panel (b) of Fig. [T] shows the 
conduction band diagram under the assumption that the 
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latter is the dominant source of dopants. 

Transport properties of the sample are measured us- 
ing ac four point method in van der Pauw geometry to 
eliminate the influence of the contacts resistance^-, and 
the film is patterned in a clover-shape. Standard lock-in 
techniques are used for the current source and voltage 
measurements. The sample is tested in an Oxford VTI 
helium gas flow cryostat, varying the temperature from 
300 K down to 1.5 K. A superconducting magnet supplies 
a magnetic field up to 15 T. 

As previously described, the temperature range is di- 
vided into three regimes according to the behavior of 
the high mobility layer in the highest magnetic fields, 
namely quantum Hall (QH), Shubnikov de Haas (SdH) 
and Drude (D) regimes. The conductance parameters of 
the high mobility quantum confined layer will be labelled 
Q, and those of the parallel layer with the symbol ||. 

We begin sample characterization at base temperature 
in the QH regime and follow the notation of Grayson and 
Fischer— for modeling parallel conduction in this regime. 
At 1.6 K, magnetic field B dependent measurement of 
the total longitudinal sheet resistance p^^J: and transverse 
Hall resistance p*^y shows oscillations in p*^* which do 
not quite reach zero in Fig. [D^c), and plateaus in p*°y 
which likewise do not quite reach the nominal quantized 
plateau values in Fig. HJa). Both behaviors are indica- 
tions of possible parallel conduction. In the low B region 
a single series of SdH oscillations periodic in 1/B begins 
around 0.5 T with no evidence of beating in Fig. [2][b), 
revealing that only one type of high-mobility carrier ex- 
ists in this system, and eliminating the possibility that 
a quantized second subband in the quantum well is the 
source of any parallel conduction^". Thus we can safely 
assume that 1.6 K is in the parallel conducting QH regime 
and follow the analysis of Ref.— . The carrier concentra- 
tion of the quantum confined 2DEG n'^ can be calculated 
from quantum oscillations to be n'^ = 2.1 x lO^^cm"^, 
and the mobility can be calculated assuming the conduc- 
tivity of the parallel channel is negligible compared to 
that of the 2DEG. The result reveals a high mobility of 
160, 000 cm^/Vs. According to Ref.-, the small carrier 
density and mobility in the parallel channel is 
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FIG. 2: (a) Plot of measured pty* and p^°' at 1.6 K. Quan- 
tized Hall conductance values h/ve^ are indicated with short 
lines and indexed by the filling factor v. Carriers are n-type 
as revealed by the sign of Hall voltage, (b) Shubnikov de 
Haas oscillations in low magnetic field, (c) p*°' = 84 $7 at the 
minimum u = 1. 



transverse conductivities ct^^ = and a'^y = ne^ /h, re- 
spectively. At higher temperatures tr^^ ^ 0, and the high 
mobility layer enters the SdH regime wherein no analy- 
sis can separate the the conductances of the two layers. 
When (T^^ is finite, the parallel channel transverse resis- 
tance pUy becomes 
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where is the total resistivity tensor. is the ma- 



trix determinant, |p*°*| = det(p*°*). By Taylor expand- 
ing Eq.Q for small cr^^. we can get a sufficient condition 
for the small longitudinal sheet conductivity limit for the 
high mobility layer where Eq.([T]) is valid: 
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where and are the measured longitudinal and 
transverse sheet resistivity, v the filling factor, e the elec- 
tron charge and h the Planck constant. These equations 
result in nil = 5.9 x 10^ cm^^ and /^H = 140cm^/Vs. 
The two orders of magnitude lower density confirms the 
assumption that parallel channel conduction is compar- 
atively small. 

Equations [1] and [5] are only valid when the high mo- 
bility layer is in the QH regime, with longitudinal and 



We can now define the SdH regime for characterizing 
parallel conduction. In the data of Fig. 3, both p^°^ and 
p*°y are shown for all temperatures below 17 K where 
oscillatory behavior is observed. For v = 1, Eq. @ is 
satisfied only at 1.6 K, confirming that this data is in the 
QH regime, and all higher temperatures from 3.8 K to 17 
K are in the SdH regime. 

Above 20 K there is no oscillatory behavior in the p^'jj 
traces of Fig. 21 however the nonlinear behavior reveals 
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FIG. 3: Measured p^"' between 1.6 K and 17 K. minima 
increase as temperature increases. Between 1.6 K and 20 K 
is the Shubnikov de Hass (SdH) regime as discussed in text. 
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FIG. 4: Magnetic field dependence of Hall resistance at tem- 
peratures above 20 K in the classical regime. Non-linear curve 
fitting is adopted to extract parameters for two parallel con- 
ducting layers. Experimental data is displayed in dots and 
fitted results in solid line. The inset shows the possible acti- 
vation energy (A_d) for the parallel conducting channel at high 
temperature above 75 K. The error bars for A_d is ±2meV. 



parallel conduction. A parallel conduction Drude model'^ 
can be applied to extract the conduction parameters for 
both the high mobility layer and the parallel channel. 
Separate sheet resistivity tensors for two layers are de- 
fined as 
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where 
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and n*,/i* and ctq are the sheet carrier density, mobility 
and sheet conductivity at zero magnetic field for layer 



i, respectively. The corresponding conductivity tensor 
obtained by inverting Eq.Q is 
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The total conductivity tensor of the system is 
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Inverting cr*°* gives the total resistivity tensor p'°* = 
(cr*°*)~^, whose components p^°^ and p^°y are directly 
measured in a four-terminal resistance or Hall measure- 
ment. p*'°l^ p*'°y are functions of B with four parameters 

cr^, (jg, [fi and /^ll. In particular when i? — OT, 
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These are the well-established expressions for mixed con- 
duction in low magnetic fields^^. With constraints given 
in Eq. ([9]), (fTOl) . a standard curve fitting program such 
as Auto2Fi1^ is used. Then carrier density is calculated 

Figure [Sj^a) summarizes the temperature dependence 
of Hall mobility and Fig. Efb) carrier density in the en- 
tire temperature range from 1.6 K to 300 K. Because 
below 20 K the parallel channel almost freezes out, rfi 
and p^ are estimated from the sheet resistivity and Hall 
measurement results in the zero B limit using a single 
conducting layer approximation. The values of carrier 
density and Hall mobility from the two models are con- 
sistent and reveal reasonable trends for each parameter 
as temperature changes, p^ is almost constant of order 
of lO^cm^/Vs below 40 K, and decreases rapidly with 
increasing temperature above 40 K to circa lO^cm^/Vs 
at 300 K. /ill is temperature independent from 40 K to 
300 K. rfi is almost constant around 2 x lO^^cm"^, only 
decreasing a little from 200 K to 300 K. And n" increases 
very drastically as temperature increases. 

The triangles in Fig. [5] represent the effective mobil- 
ity and carrier density estimated from a naive single- 
layer conductivity model. This model is typically ap- 
plied in the literature to estimate sample mobility, and 
it attributes all conduction to a single layer. With this 
simple model, the room temperature 2D mobility under- 
estimates the actual 2D mobility by a factor of 2.6 and 
carrier density underestimated by a factor of 1.6. As 
temperature decreases below 140 K, however, there is no 
appreciable difference between the single layer approxi- 
mation and the estimates of the high mobility layer from 
the parallel conduction model. 

The data reported here testify to the high growth qual- 
ity of the samples. The temperature dependence of 2D 
mobility has been published by Matsuoka. et. al}^ for 
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FIG. 5: (a) Temperature dependence of Hall mobility, (b) 
Temperature dependence of carrier density. A quantum Hall 
(QH) parallel conduction model is used at 1.6 K when 2DEG 
satisfies the QH limit in strong magnetic field (QH regime) 
and a Drude (D) two-layer parallel conduction model is used 
above 20 K (D regime). The temperature in between is the 
SdH regime where the parallel conducting layer is impossible 
to characterize. The dash line marks the temperature bound- 
ary of three regimes. Parameters with superscript eff are the 
values from the traditional single conducting layer model. 



InGaAs/InAlAs heterostructures with a saturation mo- 
bility about 6 X 10'' cm^/Vs and a knee at 77 K. By com- 
parison, the mobility here saturates at 1.6 x lO^cm^/Vs 
and shows a knee 40 K. The knee in /j,'^ is known to be 
be caused by polar LO phonon scattering, and the sat- 
uration at low temperatures is interpreted in terms of 
alloy disorder scattering^"*. The higher saturation mobil- 
ity and lower knee temperature indicate that our sample 
has less alloy disorder scattering and better quality. 

The analysis shown here also allows other energy scales 
in the band diagram to be deduced, such as the dopant 
binding energy, and the heterojunction conduction band 
offset. The activation energy of the parallel conduction 
states can be estimated by fitting an Arrehnius law = 
no exp(— Au/fesT) above 75 K as shown in the inset of 
Fig. 21 where A^i is the thermal activation energy, ks the 
Boltzmann constant, no represents the number density 
of impurities. Average A^i is about 15 meV and no is 



around lO^^cm"^. 

The conduction band offset (AEc) at the base tem- 
perature 1.6 K can be estimated with the parameters 
extracted above. Since the assumed background dopant 
density of InP bulk (~ lO^^cm"^) is too low to offer 
enough electrons for 2DEG, it is reasonable to assume 
that InP layer is entirely depleted and the main source 
of donor is interfacial states between InP epilayer and 
semi-insulating InP substrate. This assumption leads to 
a conduction band structure as shown in panel (b) of 
Fig. [T] It is clear that 
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where d = 200 nm is the thickness of InP bulk, £b the 
electrical field strength in InP, A^ the binding energy 
of interfacial states, and Ao the activation energy from 
the bottom of the quantum well to the first subband. 
According to Maxwell's equations. 
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where = 12.4 and ef^y = 14 are the relative static 
permittivity for InP and lno.532Gao.467As respectiveljsi^, 
eo the electrical constant and £w the electrical field in 
the quantum well. Then with n'^ = 2.1 x 10^* cm~^ the 
conduction band decrease in InP bulk can be calculated: 
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The interface dopant density is no — 10^^ cm^^ and Ajj 
is 15 meV as mentioned above. Ep can also be calculated 
from n'5 as 



Ef 



12.3 meV, 



(14) 



where the effective electron mass in 2DEG is approxi- 
mated by that in the bulk of Ino.532Gao.467As, me* w 
0.041mo. Finally we can estimate Ao with the model 
developed by Stern^^ 
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With Ao w 38 meV, the conduction band offset can be 
estimated as AEc ~ 677 meV, about 140 meV larger 
the result obtained by Guldner et. al. in similar 
heterostructurei^. In addition to possible differences in 
alloy content, the difference in AEc niay also be caused 
by either different interface defects with differing binding 
energy, and/or different background bulk dopant density 
which would add minor curvature to the conduction band 
not included in the above approximation. 

In summary, we have presented the magnetotransport 
results for the high quality InP/InGaAs heterostructure 
and characterized them with parallel conduction models 
in different temperature ranges. The mobility and car- 
rier density for 2DEG in the entire temperature range 
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between 1.6 K and 300 K are extracted. So are those 
parameters for the parallel channel in the classical and 
quantum Hall regime. The main source of donors is be- 
lieved to be the interfacial states between InP and sub- 



strate. The methods used in this paper can be applied 
to the general characterization of heterostructures with 
parallel conduction of layers with different mobilities and 
densities. 
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